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Abstract 

We study the spectral properties of the process of explosive percola- 
tion. In particular, we explore how the maximum eigenvalue of the adja- 
cency matrix of a network which governs the spreading efficiency evolves as 
the density of connection increases. Interestingly, for networks with con- 
nectivity that grow in an explosive way, information spreading and mass 
transport are found to be carried out inefficiently. In the conventional 
explosive percolation models that we studied, the sudden emergences of 
large-scale connectivity are found to come with relatively lowered effi- 
ciency of spreading. Nevertheless, the spreading efficiency of the explosive 
model can be increased by introducing heterogeneous structures into the 
networks. 
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In the past few years, there has been intensive research on explosive per- 
colation - the abrupt development of large-scale connectivity in networks, im- 
mediately after Achlioptas et. al. brought out their percolation model with an 
unexpected sharp transition [T] . In the model, when one increases progressively 
the number of connections between nodes in a network following the suppression 
principle 0, above some critical threshold, a giant connected cluster emerges 
suddenly. The resulting transition was originally thought to be discontinuous. 
It is shown however to be continuous for some models later [3J HJ [SJ [7] . In 
addition to the explosive percolation of a single giant connected component, the 
abrupt creation of stable multiple giant components are shown to be possible 
as well [SI [H] • More recently, strongly and weakly discontinuous transitions are 
introduced to classify how explosive a percolation transition is [10]. By now, 
the study of discontinuous percolation transition has been extended to scale- 
free networks 112) , real- world networks with community structure |13j , two- 
and higher-dimensional lattices [HI [JS1 EH H3, etc- Concurrently, other ap- 
proaches like the weighted rule [T8] and the Hamiltonian approach [H] have 
also been proposed to obtain similar transitions. Most notably, this process 
has shown to be a possible mechanism that explains the growth process of a 
couple of real networks, including the human protein homology network [20] 
and nanotube clustering [3T]. For real- world systems, such behavior can have 
vital consequences when the addition of single links jTS] [22] may drastically 
change macroscopic connectivity and hence the dynamic and function of the 
networks. For instance, in a neuronal circuit, this would mean that the growth 
of one or a few additional synaptic connections might drastically alter the in- 
formation processing function in the brain. Similarly, the establishment of a 
small number of specific social relations may significantly increases the possible 
extent of infectious diseases or rumors. Here, we study the evolution of spec- 
tral properties of the network during the process of explosive percolation. In 
particular, we explore how the maximum eigenvalue of the adjacency matrix of 
the network evolves as the density of connection increases. This is of interest 
as the maximum eigenvalue of the network governs how information or diseases 
spread [231 123] . If it increases explosively due to the sudden emergence of a 
giant component, the spreading process will be enhanced. On the other hand, if 
the maximum eigenvalue remains small at the percolation threshold, an efficient 
spreading is not possible despite the existence of a giant connected cluster in 
the network. 

The smallest cluster (SC) model is one of the simplest model that show 
strongly discontinuous percolation transition. It starts with N = 2™ isolated 
nodes. At each step, the two smallest clusters in the network are identified and 
merged into a larger cluster through the creation of a link between them. In this 
case, N/2 clusters of size 2 are created during the first phase, i.e. the first N/2 
steps. In the second phase, N/4 links are added, each connecting two clusters of 
size 2 into a cluster of size 4. The process continues and by the end of the phase 
y = n— 1, only two components remain, each with size N/2. Then, the addition 
of the next edge connects these two components, resulting in a jump in value of 
N/2 in the size of the largest component (S). If the nodes to be connected in 
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each step are chosen randomly from the two smallest clusters, then at the end 
of phase y, since N/2 V edges are added to the network and the total degree is 
increased by N/2^ y ~ 1 \ the evolution of the degree distribution will be described 
by: 

8P(k, y)/dy = ^zi{-P{k, y) + P(k - 1, y)) . (1) 
With this, we can obtain the probability of nodes with degree y at phase y: 

P{y,y) = ^r-.. (2) 

Note that P(k,y) = for all k > y. In the earlier phases, the maximum degree 
of the network, k m increases linearly with the phase. Nonetheless, since P(y, y) 
drops rapidly as y increases, beyond a critical phase, the probability for k m to 
increase becomes very small. Here, we define y c to be the critical phase where 
P(y c ,y c ) — 1/N. The value of y c can then be obtain by solving the following 
equation: 

y c (y c - 1) = 2n. (3) 

Given that the maximum eigenvalue (A m ) of a network's adjacency matrix is 
inversely proportional to the square root of the network's maximum degree, A m 
is not expected to grow explosively at the percolation threshold as the size of 
giant component does. In fact, as shown in Fig. 1(a) for a smaller network 
size and Fig. 1(b) for a larger network size, the maximum eigenvalues do not 
increase further after the critical phases. Note that the results are obtained 
through numerical iteration of Eq. ([1]) with the maximum eigenvalue being 
approximated as the maximum eigenvalue of the largest cluster |25) : 



A "=f- + !P < 4 » 

For the smaller network size, a result based on an averaging of 20 network 
realizations is also obtained for comparison. 

On the other hand, if the nodes to be connected in each step are chosen 
to be the largest-degree nodes in the two smallest clusters, then at the end of 
phase y, 

forfc = fe m + l, (5) 
0, for all other k . 

In this modified SC model, the maximum degree increases linearly with the 
phase. Thus, as shown in Fig. 1, larger maximum eigenvalues and more efficient 
spreading are obtainable at the percolation threshold for giant clusters that 
emerge explosively. 

Next, we look at a more general model, the Gaussian model In this 

model, the distribution of the cluster size is controlled to be Gaussian by im- 
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plementing the following probability of occupation: 



p = exp 



(6) 



for an edge that is selected randomly among the empty ones. In the equation, s 
denotes the size of the cluster that would be formed by adding the selected edge 
while s denotes the average cluster size after adding the edge. When an intra- 
cluster link is chosen, s will be taken as twice the cluster size. The parameter 
a controls the size dispersion and we consider only for cases where a > 0. For 
a = 0, all edges have the same probability to be added and the model gives the 
Erdos-Renyi graph. In this case, a cluster of relative size is expected to appear 
at the percolation threshold t c = 0.5 and grows gradually as the connection 
density increases (see the inset of Fig. 2). For a > 0, the formation of a cluster 
with size which differs significantly from the average cluster size is suppressed. 
The result is a sudden emergence of a giant connected cluster at t c ~ 1. In the 
inset of Fig. 2, we show the explosive percolation transitions for a = 0.5 and 
a = 2. 

For this model with a > 0, we consider only the addition of inter-cluster con- 
nections for the derivation of the evolution of the network distribution. Letting r 
to be the ratio of the total number of sampled edges to the network size, then at 
t, the sampling probability for an edge between a node in the cluster of size s and 
a node in the cluster of size s' is 2Q(s, t)Q(s', t) where Q(s, r) = ^ fe P(s, k, r). 

If this edge is 



The probability for this edge to be added is exp 

added, P(s,k,r) will decrease by sP(s, k,T)/NQ(s,r). On the other hand, to 
form a cluster of size s, the selected nodes must be in clusters of sizes u and v 
with u + v = s. The probability for a connection to be added between these 
two nodes is exp [—a ( s -= s -)] ■ If this edge is added, nodes with degree k in the 
two clusters will increase the fraction of nodes with cluster size s and degree k 
unless they are the selected nodes. Then again, nodes with degree k — 1 in the 
two clusters will increase the fraction of nodes with cluster size s and degree k 
if they are the selected nodes. The evolution of the distribution of the network 
can thus be written as: 

- -2^M,^Q(/,r)eh(^)] 

s' 

+ £ Q(u,T)Q(v,T)e[-<^)] 

u-\-v—s 

(u-l)P{u,k,r) P(u,k-l,r) 
NQ(u,t) + NQ(u,t) 

(v-l)P(y,k,T) P(v,k-1,T) \ 

NQ(v,t) + NQ(v,t) J " 1 ' 

Note that the ratio of the total number of added edges and the network size can 
be approximated as t — (fc)/2. 
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Equation ([7]) is iterated numerically, then the degree distribution R(k, r) = 
J2 S P( s > k, t) is used to calculate X m for r with J^k T ) = !• The growth of 
X m versus t is shown in Fig. 2 for both the random network and the explosive 
model. As shown, A TO increases gradually as t increases. Interestingly, while the 
size of the largest connected component at t 1 is considerably larger for the 
explosive models, the maximum eigenvalue is found to be consistently smaller 
throughout the process. In other words, although the connected component 
is larger in networks generated through the explosive model, the efficiency of 
spreading in these networks is lower compared to that of the random graphs. 
In addition, the sudden emergence of a large-scale connectivity in the explosive 
model comes at a cost of lowered navigation efficiency in the network. In Fig. 3, 
the growth of average inverse path length (L) are shown for networks developed 
through the Gaussian model with a = 0, 0.5 and 2. For a > 0, although the 
average inverse path length increases explosively at the percolation threshold, 
it is still significantly smaller than that of the random networks. 

In order to increase the spreading efficiency of the networks, we implemented 
a modified Gaussian model. At each time step, instead of drawing two nodes 
randomly from the network, nodes i and j are drawn from the network with 
probabilities w,;/ J2m Wm an< ^ w i I J2 m Wm respectively. Here, w m are weights 
assigned to nodes in the network following the Chung and Lu (CL) model [TT1|2"6] 
which is used to construct artificial scale-free networks. Specifically, the nodes 
are indexed as i — 1, 2, • • • ,N and each node is assigned a weight of 

(i + io-1)-" (8) 

for fj, < 0.5, 
for0.5 <fi< 1. 

Then, if edge ij does not already exist, it will be added with the probability 
described by Eq. ([5]). Numerical simulations are performed for this modified 
model. Indeed, as shown in Fig. 4, X m is greatly increased when heterogeneous 
structures are introduced into the explosive model. 

Explosive percolation model has altered the traditional understanding on 
percolation transitions. While the sudden emergence of a large-scale connectiv- 
ity may significantly change the dynamic and function of a network, for infor- 
mation and mass transport in real-world networks, the existence of a large-scale 
connectivity alone is not sufficient. The fact that most of the real- world net- 
works are small- world networks, in which the average path length is scaled as 
log(-ZV), reveals that navigation efficiency in a network is a critical consideration. 
However, as shown by our results, the sudden emergence of a giant component 
in the explosive model comes at a cost of lowered spreading efficiency in the 
network. As such, the insertion of heterogeneous structures into the network 
becomes important in achieving high spreading efficiency in explosive models. 
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Figure 1: The maximum eigenvalue at the end of each phase for the smallest 
cluster model (lower curves) and the modified smallest cluster model (upper 
curves) of size (a) N = 2 13 and (b) N = 2 55 . Note that the ensemble averages 
obtained by simulation of the SC model with 2 13 nodes are shown as triangles 
whereas the approximations obtained through numerical iteration of Eq. ([1]) 
and ([5]) are shown as circles. 
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Figure 2: The maximum eigenvalue versus the ratio of the added edges for the 
Gaussian model with a = 0, 0.5 and 2. Note that the numerical results are 
obtained through an average of 20 network realizations of size N = 2 13 and 
are shown as circles (a = 0), squares (a = 0.5) and crosses (a = 2) whereas 
the analytical approximations are shown as dashed (a = 0), thick (a = 0.5), 
and solid (a = 2) lines. For numerical iteration of Eq. (O, the parameter used 
are s = 1, 2, ■ ■ • ,N, k = 0, 1, ■ • • ,25 and < t < 0.85. Insets: The percolation 
transition for the Gaussian model with a — (solid line), 0.5 (line with squares) 
and 2 (crosses). 
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Figure 3: The average inverse path length versus the connection density for 
the Gaussian model with a = (thick line), 0.5 (squares with dashed line) 
and 2 (crosses). Note that results are obtained by averaging over 20 network 
realizations with TV = 2 13 . 




Figure 4: The maximum eigenvalue versus the ratio of the added edges for the 
Gaussian and modified Gaussian model. Note that the numerical results are 
obtained through an average of 20 network realizations of size N — 2 13 and are 
shown as squares (Gaussian model with a — 0.5), triangles (modified Gaussian 
model with a = 0.5), crosses (Gaussian model with a = 2) and plus signs 
(modified Gaussian model with a = 2). 
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